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Abstract 

A weak Guderley-Morawetz problem is formulated for a mixed 
elliptic-hyperbolic system that arises in models of wave propagation 
in cold plasma. Weak solutions are shown to exist in a weighted 
Hilbert space. This result extends work by Yamamoto. MSC2000: 
35M10. 



1 Introduction 

A characteristic feature of wave propagation in cold plasma is the possibility 
that a hybrid resonance surface, along which the linearized equation for the 
electric field changes from elliptic to hyperbolic type, may be parallel to a flux 
surface. This property can be represented in two dimensions by setting the 
hybrid resonance curve tangent to the line x = at the origin of coordinates. 
The situation is somewhat different from that found in, for example, linear 
models of transonic fluid dynamics [eq. (3), below]. In that case the sonic 
line is everywhere normal to the line x = 0. 

A model for such a resonance curve is the equation 

x = a{y), 
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where a (y) is a continuously differentiable function of its argument satisfying 

a (0) = a' (0) = 0. 

In addition, we assume for simplicity that both a {y) and a' {y) exceed zero 
for y exceeding zero. 

This leads us to consider mixed elliptic-hyperbohc systems having the 
form 

Lu = /, (1) 

where 

u = {ui (x, y) , U2 (x, y)) , / = (/i (x, y) , {x, y)) , 
{x,y) en gRx {R/R-) , 
(Lu)^ = [x - (7 (y)] ui^ + Kui + U2y, 

(Lu)^ = Uiy - U2x 

with data 

dx dy „ 
ui— + U2-f-^0 2 

as as 

given on a portion of the boundary of fl. Here K is a. constant in [0, 1] and 
ds denotes the line element on dQ. The system is eUiptic for x > cr (y) and 
hyperbohc for x < cr (y). Following [MSW], we emphasize the analogy to 
fluid dynamics by calling the curve x = a (y) the sonic curve. 

In the cold plasma literature, eqs. (1), (2) tend to appear in scalar- valued 
special cases. In all these cases a (y) is proportional to y^, but this specific 
restriction is not imposed by the physics; concerning the physical model, see 
[Wl], [W2]. U ui ^ tp^, U2 ^ ipy, a- (y) = y^ and / = (0, 0), the system 
reduces to a scalar equation introduced in Sec. 3 of [MSW]. In the context 
of this equation, condition (2) corresponds to imposing constant boundary 
conditions on the scalar solution ip {x, y). A uniqueness theorem was proven 
in [MSW] for X = 1/2, in order to show the existence of a domain on which 
the classical Dirichlet problem is ill-posed for the equation. Numerical argu- 
ments for a complex perturbation were also introduced. A similar equation 
[cT (y) oc y^, Ui = ip^, U2 = —'ipy, K = 1, f = (0,0)] appeared earlier in the 
physics literature, also in the context of wave propagation in cold plasma 
[PF]. In this case certain exact solutions were constructed. Finally, sys- 
tem (1), (2) in the case a (y) — y^ was studied in an interesting unpublished 
Ph.D. dissertation [Y] on the existence of weak solutions possessing Dirichlet 
data on a certain small domain near the origin of . 
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These equations share with the Tricomi equation 

yi'xx +i'yy = ^ (3) 

a multipUcity of possible approaches to formulating boundary-value prob- 
lems. From a mathematical point of view, the Dirichlet problem, in which 
data are assigned on the entire boundary, is the "wrong problem" to solve for 
equations of mixed type, as this problem tends to become over-determined 
in the hyperbolic region. Well-poscd problems for elliptic-hyperbolic equa- 
tions generally include a characteristic gap on which data have not been pre- 
scribed. In addition to the example of an over-determined Dirichlet problem 
for the scalar equation considered in [MSW] , there are analogous examples 
for Tricomi-like equations [M2]. 

However, physical applications of eq. (3) to transonic fluid dynamics, 
and of scalar forms of (1) to wave propagation in cold plasma, suggest that 
it should be possible to prescribe data over the entire boundary. This con- 
tradiction suggests that classical solutions will have httle apphcation to such 
physical problems. In terms of weak solutions u to the system (1), the 
Dirichlet problem requires solutions to lie in an appropriate Hilbert space Hi 
and satisfy 

- (m, LV)^2 = (/, 

V</? = i'^iT^'i) G where H2 is another appropriate Hilbert space and 
where (/^^^ = on the entire boundary of the domain [c.f. [M4], eqs. (6)-(8), 
for the Tricomi case]. This problem may or may not be well-posed. A 
different approach is to require that the component of the test function (/? 
vanish only on the noncharacteristic part of the boundary, and that </? satisfy 
condition (2) on characteristics. In the early literature (see, e.g., Sec. 4 
of [M3]) an elliptic-hyperbolic problem having Dirichlet data given on the 
entire boundary is called the closed, or full Dirichlet problem to distinguish 
it from the mathematically natural case of Dirichlet data given off of the 
characteristics. Following [LP], we prefer instead to distinguish the prob- 
lem in which data are given on the noncharacteristic part of the boundary 
by calling such a problem Guderley-Morawetz, reserving the term Dirichlet 
problem for Dirichlet data given on the entire boundary as in [M4] and [P]. 

The existence of weak solutions to a Guderley-Morawetz problem is proven 
for a Tricomi-hke system in [Ml]. The estimates in Sec. 2.5 of [Y] can be 
extended to imply the existence of weak solutions to a Guderley-Morawetz 
problem for eqs. (1), (2) on a relatively large and general domain. This is 
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the content of Theorems 2 and 3 of Sec. 3, below. The arguments in [Y] 
assume that the weak Dirichlet problem and the weak Guderley-Morawetz 
problems are identical. (Concerning this identification, see the Remark fol- 
lowing the proof of Proposition 1, below.) It is claimed there, on the basis 
of Guderley-Morawetz estimates modeled on [Ml], that weak solutions of 
(1) exist for Dirichlet data prescribed on the entire boundary. We make no 
such claim for the generalization of those estimates given here. However, 
the techniques used to prove the existence of weak solutions to a Guderley- 
Morawetz problem for (1), (2) will yield a uniqueness theorem for strong 
solutions to the Dirichlet problem for this system, on a more restricted do- 
main, almost "for free" (see [Y], Theorem 1, Sec. 2.7; also see our Remark, 
Sec. 3.3). Moreover, it is possible to derive the existence of weak solutions 
to the Dirichlet problem for eq. (3) by considering a sequence of Guderley- 
Morawetz problems in which the characteristic gap is "marched" to a singular 
point on the sonic curve; see [G], [M3], and [P]. It is a reasonable conjecture 
that this method can be modified to apply to systems such as (1), (2) on an 
appropriate domain, but this is not attempted here. 

Equations (1) cannot be mapped into a system of the form studied in [Ml] 
on any domain that includes the origin. On the one hand, relatively little is 
known about such elliptic-hyperbolic systems which do not directly general- 
ize the Tricomi equation. On the other hand, the method of proof adopted 
here is by now quite standard. It is required to find a Hilbert space U, a 
domain Q, and a multiplier M under which weak solutions can be shown to 
exist, without unreasonable restrictions on generality, by a uniqueness-plus- 
projection argument using the abc method and the Riesz Representation 
Theorem. Because this system comes from a physical model, we addition- 
ally hope that our conditions on U and Q will be physically reasonable. 
For instance, physical/numerical arguments for special cases suggest that a 
singularity should be permitted at the origin [MSW]. This influences the 
weighting of the Hilbert space U, as does the existence of particular physical 
solutions (Sec. 2.2.2). 

We note that every nondegenerate conic section is equivalent under the 
projective group to the unit circle. In that sense, the system (1), (2) with the 
choice a (y) = is "gauge equivalent" to a system on the extended projective 
disc studied in [O] by similar methods. In that system, the elliptic part of 
the domain has a geometric interpretation as hyperbolic points in projective 
space and the hyperbolic part, as ideal points. As it is not clear that 
projective invariance has any physical meaning in the context of cold plasma 



4 



dynamics, this analogy will not be pursued. 



2 Formulation of the boundary-value prob- 
lem 

2.1 Domain 

In proving weak existence for the Guderley-Morawetz problem we assume 
that the domain Q, having piecewise continuous boundary, is enclosed by the 
arbitrarily large but finite rectangle 

R = {{x,y) \ — oo < p < X < £, < y < q < oo} , 

where i, p, q are fixed but arbitrary real constants. We assume that R has 
been chosen so that the distance along the x-axis from sup^ x to £ is an 
arbitrary but fixed positive number S. Because we are assuming the existence 
of an elliptic region for eqs. (1), we take £ > 0. 

The elliptic region of Q consists of the region of the first quadrant bounded 
by the sonic curve x = a (y) and a smooth curve Ci emerging from the origin, 
along which 

dy 

> 

dx 

with equality only at the origin, and 

a{y)^ + b{y)<0 (4) 

for specified functions b{y) < and a{y) > 0. Wc assume that Ci intersects 
the sonic curve at a point {xQ,a{xo)) G fi, where Xq > 0. For example, if 
a, 6, and a are defined as in Theorem 3 of Sec. 3, then the family of curves 
given by y = ex'^ for m > 1/2 and x >0 satisfies condition (4) for 



n-2m 

e < 

Km 



whenever K > 0. [Condition (4) is automatically satisfied for such a, b, and 
a ii K — 0.] If we further specify m < 1/K, then we are guaranteed that 
Xq < i — S provided we choose 



l-2m 
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The hyperbolic region is bounded by a piecewise smooth curve F U C2, 

where: 

r is a characteristic curve 

^ = -V'^iy)-^ (5) 

emerging from the sonic curve at {xq, a (xq)); 

C2 is a piecewise continuous curve intersecting the characteristic F at a 
single point on the left endpoint of C2 and intersecting Ci at the origin on 
the right endpoint of C2. We assume that dy < and dx > on C2. We 
orient the boundary in the counterclockwise direction. 

The relation of this domain to the domain considered in [Y] is discussed 
at the end of Section 3.3. 

2.2 Function spaces 
2.2.1 WeEik solutions 

Denote by U the Hilbert space consisting of all pairs of measurable functions 
{ui, U2) such that 



\u\ 



cr' {y) {ul + K2) dxdy 

.J Jn 



1/2 



is finite. Here 

(u, w)^^ I I a' (y) (uiwi + U2W2) dxdy. 




Denote by W the linear space of continuously differentiable functions (wi, W2) 
vanishing at the origin of and satisfying: 

Widx + W2dy — (6) 

on the characteristic F, toi = on dVL/V, and 

1/2 



< 00, 



where 

{L*w)^ = [x - cr {y)] wi^ + {1- K)wi + W2y 
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and 



{L*w)2 = Wiy - W2x- 



We define a weak solution to eq. (1) under the boundary condition (2) to 
be any u & U such that Vw e W, 



(w, /) = - {L*w,u) 



(7) 



under the inner product ( , ). 

Denote by H the Hilbert space of measurabie functions {hi, /i2) for which 
the norm 

1/2 



is finite. The inner product on H is given by 

{h, 9)* ^ j J -^r^ {hiQi + ^252) dxdy. 
The prescribed data / will be assumed to lie in the space H. 



2.2.2 Similarity solutions 

Analysis of scalar special cases of system (1) with a {y) = y^ suggests the 
presence of a singularity at the point x — y — 0. See, for example, [MS W] , 
where this is discussed in detail. This would suggest a radial weight for the 
energy functional. Quadratic radial weights are applied in [Y]. However, 
the space U constructed here also arises naturally in connection with eqs. 

(!)• 

As a simple example, consider similarity solutions for the case a (y) = y^ 
having the form ui — ip^, U2 — ipy, and 

where 1/ is a parameter and F satisfies the hypergeometric equation 
(1 - /x) {v-l)F (/.) - 2 (z. - 1) /.F' (/.) + /.^F" (/.)] 
+ [2F' (pi) + ApiF" (/.)] = 
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for 




Properties of such solutions for complex values of v arc studied in [MS W] . 
We consider here the case of real- valued v, as in [PF] and [Wl]. It has been 
observed [Wl] that if \x\ is sufficiently small, then F jjy ox F 
Taking F ~ /x*" for i/ = 1/4, we find that 



If u lies in the function space U ^ then ij) has weighted Dirichlet norm 



In fact, Eu {i>) is finite on Q Vi/ > 1/4. If we include solutions of the form 
F ~ fj,"-'^, then Eu would be finite on Q Vi/ > 5/4. 

2.3 The weak problem is well-posed 

Proposition 1 Any twice continuously differentiable weak solution of the 
Guderley-Morawetz problem for (1), (2) on Q, as defined by eq. (7), is a 
classical solution. 

Proof. We refer to the domain as Jl, but the argument also holds without 
alteration on much more general domains. For u & U and w integration 
by parts yields 

{u,L*w)= / / ui {[x - a (y)] + - K) Wl + W2y} dxdy 








8 



- / {wiU2 + W2U1) dx+ 
Jon 

{[x-a{y)]wiUi-W2U2}dy. (8) 



I an 

On dfl/r, wi = 0, implying that 



u,L*w):g^,^ = - W2 {uidx + U2dy) . (9) 
Jan/r 



'an/r 

Equations (5) and (6) hold on F, implying that 



, I/w)|p = J - {wiU2 + W2U1) dx + {[x - a (y)] wiUi - 'W2U2} dy 
= J ui i^[x — a (y)] wi — ""^2^ I dy — U2 {widx + W2dy) 



x-a{y) + 



dx^ ^ 



wiUidy = 0. (10) 



Substituting eqs. (9) and (10) into (8) and using (7), we obtain 

— {w, /) = (m, L*w) = — {Lu, w) — W2 {uidx + U2dy) . 

Jan/r 

Because this identity must hold for every w E W, we conclude that the 
quantity Uidx + U2dy must equal zero almost everywhere on dD,/r. Applying 
the hypothesis that u is twice continuously differentiable, we complete the 
proof of Proposition 1. 

Remark. The value of the 1-form uidx + U2dy on the characteristic F 
is left undetermined by a definition of weak solution based on (7), so this 
argument will not establish the well-posedness of the weak Dirichlet problem 
for (1), (2) on fl (unless we change the boundary conditions on w to wi = 
on dQ). However, classical solutions u of either the Dirichlet problem or 
the Guderley-Morawetz problem satisfy (7). This ambiguity seems to be the 
basis for the attempt in [Y] to identify the weak forms of the two problems; 
c.f. Sec. 2.3 of [Y]. 
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3 Results 



Theorem 2 Let K e [0, 1/2]. Let the functions a{y) and b{y) in condition 
(4) be given by 

-1 



a{y) = K 



y + 



and 



b{y) = - 



1 + 



a (t) dt 



For every f E H, there exists on fl a weak solution to system (1) with the 
boundary condition (2) given on dfl/F. 



Theorem 3 The conclusion of Theorem 2 extends to the case K e [0, 1] if 
we replace the definitions of a{y) and b{y) by 



and 



and specify a (y) — y^. 



aiv) = Ky 



1 + 



2£ 



The proofs of Theorems 2 and 3 modify the argument in [Ml]. In 
addition, we adapt a number of choices made in [Y] , which is also based on 
[Ml] . The resuhs follow from an a priori estimate. 

Lemma 4 Under the hypotheses of either Theorem 2 or Theorem 3, 3k > 
3 Vw E W we have 

k \\w\\^ < \\L*w\\* . 



3.1 Proof of Lemma 4 

We prove the lemma by the abc method. Let 



M 



a b 
c d 
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where a and b are given by the hypotheses of the theorems; c and d will be 
chosen. Then 

/ = {L*w, Mw) = 

/ / {[x - a {y)\wix + {I- K)wi + W2y}{awi + hw2)dxdy 
J Jn 

j I i'^iy ~ '"^2a;) {cwi + dw2) dxdy. 

Notice that a and b are defined so that ax = bx = 0. We find as in [Y], Sec. 
2.4, that 

a[x-a (y)] WiWix =^{{a[x-a (y)] wl}^ - awf) ; 

bw2 [x- a {y)] wix = 
{b[x-a (y)] WiW2}^ - bwiW2 -b[x-a (y)] WiW2x, 

awiW2y = {awiW2)y - ^ {awf)^ - ayWiW2 + aw2W2x - awiyW2] 
bw2W2y = ^ {bwl)y-bywi\ ; 



2 
1 
2 

Choose 

d — a 



cwiyWi — - 



{o^Tjy - Cywi^ 



and c — —[x — a (y)] b. Taking into account cancellations, we can write 
I — Ii + I2, where I2 is a line integral and 

h — [awf + 2/3wiW2 + 7^2) dxdy 

J Jn 



for 



(^ = l{by[^-<^ (2/)] - b {y) a' {y)} + Q - i^^ a 



iy), 



and 



P = -^[ay + Kbiy)], 



1, 

7 = -^by. 
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Case 1: Under the hypothesis on K in Theorem 2, the coefficient of a{y) 
in a is nonnegative, and we can write 



and 



Thus in this case we have 



7 = 



y + i~^ / a{t)dt 
Jo 

5a' (y) 

/3 = 0, 



+K{--K 



> 



2e 



Ii> J j {awf + 7^2) dxdy > j j iv) (^'i + tJ^l) d^^V^ 

where x = {^A}- 

Case 2: Under the hypotheses of Theorem 3, we have 



a 



y_ 
2e 



{2y^ + 2i-x) + (^^-KjKy 



1(2,^ + 2.-.)-!^ 
y-(2y^ + i-x)>X 

2i^ ^ ' - 2e 



> 



/3 = 



and 

Notice that 



^ 2t 



2 



.2iJ 



-x + -y 
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Cauchy's inequality implies that 

2(3wiW2 > —2 \P\ \wi\ \w2\ > —1\fa \w\\ \w2\ > —awl ~ I'^l (12) 

in fi/ = 0} . This already implies that the W^-norm of w is positive inside 
the upper half-plane. It remains, however, to derive an explicit lower bound 
on the coefficient of y{w\ + 

We claim that there is a constant e e (0, 1) , depending only on R, for 
which 

0<a7-5(^)'<£a7. (13) 

To estabhsh this claim, note that the left-hand inequality in (13) is obvious 
from (11), and the right-hand inequality will be satisfied provided 

«7(l-e)<i(J^;. 

Assuming without loss of generality that y exceeds zero (the inequality is 
true trivially otherwise), our criterion becomes 



2 + £ 



l+[]^-K\K 



— X < 



Replace the quantity on the left by its largest possible value, given that 
(1/2 — K)K < 1/16. Our requirement becomes that e be chosen sufficiently 
close to 1 so that 

2q^ + — -p< 



l-£' 

or 

The quantity on the left exceeds zero, as £ — p exceeds 6. 

Now (11) and (13) imply that inequality (12) can be improved to read 



2PwiW2 > -2 \P\ \wi\ \w2\ > -2^ aj-S |wi| 1^21 > 



—2y/ecrf \wi\ \w2\ > —^a.w\ — \fe^w. 

Thus in this case 



/i > (l — \fe) I I [awl + 7^2) dxdy > 
J Jn 



13 



''y {Swl + wi) dxdy > / I y {wl + wl) dxdy. 




2£ J Jn ' ' 2£ 
The remainder of the proof is identical for either set of hypotheses. The 
boundary terms resulting from applying Green's Theorem on f2 are given by 

- J ^- {wl - [x - a {y)] wl} + awiW2^ dx+ 

/ ( ^ { - ^ (y)] wl -wj} +b[x - a (y)] W1W2) dy. 
On C2, = 0, > and dy < 0, so the signs of a and h imply that 



h\C2 ^-T, j aw'idy + bw'^dx > 0. 
^ -'C2 



On the characteristic T we have x < a (y) , dx < 0, and % > 0. In 
addition, eqs. (5) and (6) imply that 

wl + [x-a (y)] wj = wf [a (y) - x] + [x - a (y)] wf = 0. 

We have 



where 



and 



h\r — hi + -^22, 
/ awiW2dx + - f a {[x — a (y)] wl — wl} dy 
^{wl + [x-a {y)] wl} dy^O 

-^22 = - y ^{wl-[x-a (y)] wl} dx 
+ J b[x-a (y)] wiW2dy = 

-{wl + [x-a{y)] wl}dx^O. 

On Ci inequality (4) holds; in addition, dx > 0, dy > 0, and Wi = 0. 
Writing 

.dy 



aw2dy = / aw2--j-dx, 
Ci Jci dx 

14 



we find that 



l2\Ci I awjdy + bw2dx = 



Ci 

Wndx. 



a— + 

ax 



'Ci L 

Inequality (4) implies that the integral on the right is nonnegative. 

The preceding arguments establish the lower bound of the lemma. 

In order to obtain the upper bound for the inequality of Lemma 4 we 
reason in either Case 1 or Case 2 as in [Ml] , writing 



lim f , ^ ^ {^^Ja' (y) + rj Mv}j < 



L*w 



C{M)\\L*w\\*\\w\\^. (14) 

The constant C{M) will be a finite positive number provided the functions 
a, 6, c, d are bounded. The existence of such a bound follows from the finite 

character of the constants p, g, and C. 

We obtain, under the hypotheses of either theorem, the inequality 

C'{p,q,l,e,5)\\w\l<\\L*w\\* 

for C > 0. This completes the proof of Lemma 4. 

3.2 Proof of Theorems 2 and 3 

Both theorems follow from Lemma 4 by a standard argument. An inequality 
similar to (14) implies that Vw e VF, 



\{wJ)\<c^\\L*w\\ 
for a constant Cq depending only on For fixed f & H, the functional 

G{L*w) = {w,f) 

on L*w can be extended to a bounded hnear functional on H. The Riesz Rep- 
resentation Theorem then implies the existence of an element h — {hi, /i2) £ 
H for which 

{wj)^{L*w,hr. 
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Defining u — {ui, U2) , where 



and 



U2 = 



we find that u E U and 



iwj) = {L*w,hr 



{L*w, u) 



Vw e which completes the proof. 
3.3 Remark 

By slightly modifying the proof of Lemma 4 it is possible to prove the unique- 
ness of strong solutions to a Dirichlet problem on a more restricted domain. 
Replace ^2 by a domain in which C2 is replaced by the piecewise linear 
curve Ai U A2, where Ai is a vertical line segment x — const. < 0, lying in 
the interior of R, bounded above by F and below by A2; such vertical lines 
correspond to flux surfaces in the cold plasma model; A2 is the segment of 
the X-axis bounded on the left by the line segment Ai and on the right by 
the line x = 0. The curves Ci and F are identically defined on Q and fl'. 
Let condition (4) be satisfied for b{y) defined as in Theorem 2 and for 



Let K lie in the interval [1/2, 1] . Then for every f E H there exists at most 
one strong solution in U to eqs. (1) on Q.' with the boundary condition (2) 
given on almost all of dO.' . This conclusion extends to the case K G [0, 1] if 
h{y) and a (y) are defined as in Theorem 3 and a (y) = {1 — K) y. 

By a strong solution of (1) we mean an element u & U ior which there 
exists a sequence u" & U such that 




lim Wu'' - 



u 



* 







and 



lim WLu" - /II* = 0. 
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This strong solution satisfies the boundary condition (2) on almost all of dQ' 
if in addition 



where ds is the line element on dO,'. 

Suppose that we impose the following additional restrictions and modifi- 
cations on the domain Q' : 

the arbitrarily large finite rectangle R in the upper half-plane is replaced 
by a sufficiently small circle Rq in the upper half-plane, tangent to the origin; 

the hue segment Ai is chosen to he sufficiently close to the y-axis; 

the line segment A2 is replaced by that segment of Rq bounded on the 
left by Ai and on the right by the y-axis; 

the characteristic curve F is a curve satisfying 



emerging from the parabola x = aX a. point [5,5 j sufficiently close to the 
origin; 

condition (4) is satisfied for h{y) = —(1 + y"^) and a{y) = (1 — K)y. 

Then Q,' becomes identical to the domain D considered in Ch. 2 of [Y]. 
A uniqueness theorem for solutions of (1), (2), with a [y) = y'^, lying in a 
radially weighted Hilbert space over D is given in Sec. 2.7 of [Y]. 
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